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Abstract. We provide a model to investigate feedback control of entanglement. It consists of two distant 
(two-level) atoms which interact through a radiation field and becomes entangled. We then show the 
possibility to stabilize such entanglement against atomic decay by means of a feedback action. 

PACS. 03.67.Mn Entanglement manipulation - 42.50.Lc Quantum fluctuations 
1 Introduction idea of Ref.jlj to a recently proposed model |S] consisting 



Over the last decade, the rapid development of quantum 



of two distant (two-level) atoms (two qubit) which interact 
through a radiation field and becomes entangled. We then 



_ ^ ' technology has led to the possibility of continuously mon 
K> , show the possibility to stabilize such entanglement against 

; I itoring an individual quantum system with very low noise 

atomic decay by means of a feedback action. 

and manipulating it on its typical evolution time scale Ti . 



2 The Model 



It is therefore natural to consider the possibility of con- 
trolling individual quantum systems in real time by using 
feedback. A theory of quantum-limited feedback has been We consider a very simple model consisting of two two- 
introduced by Wiseman and Milburn [21131 ■ Among recent level atoms, 1 and 2, placed in distant cavities and inter- 
developments we mention the feedback stabilization of the acting through a radiation field in a dispersive way. The 
state of a two level atom (single qubit) against amplitude two cavities are arranged in a cascade-like configuration 
damping such that, given a coherent input field with amplitude A 
Because of the relevant role played by entanglement in one of them, the output of each cavity enters the other 
in quantum processes, it would be straightforward to also as depicted in Fig^ Then, it is possible to show [SI, after 
consider its feedback control. Here we extend the basic eliminating the radiation fields, that the effective interac- 
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Fig. 1. Schematic description of the considered set-up. Two 
distinct cavites, each containing a two-level atom (1 and 2 re- 
spectively), are connected via radiation fields (solid lines). A 
coherent input of amplitude A is provided in one of them. 
Furthermore, LI, L2 and Dl, D2 represent local operations, 
namely driving fields and homodyne detection respectively. I 
is the current arising from local measurements and Fl, F2 
indicate the consequent local feedback actions (dashed lines). 

tion Hamiltonian for the internal degrees of the two atoms 
becomes of Ising type, namely 



Hint = ■iJCJi <J2 



(1) 



where crj^j^'^' indicate the usual Pauli operators. Hereafter 
we shall also use (Ti,2 = {cr[^2 + ''-'^i^l)/'^- The spin-spin 
coupling constant J scales as radiation pressure |Ap and 
goes to zero for negligible cavity detuning 

To get entanglement in an Ising model, it is necessary 
to break its symmetry 6 . To this end, we consider local 
laser fields applied to each atom (LI and L2 of Fignj such 
that a local Hamiltonian Hdrive given by 



{a e R) 



(2) 



acts in addition to Hint- Thus, the total Hamiltonian of 
the system results 



(3) 



Let us introduce the ground and excite atomic states 
|ff)i,2, |e)i.2 as eigenvectors of (j\ 2 with —1 and -f 1 eigen- 
values respectively, and rj = a/ J. Then, the eigenvectors 
of the Hamiltonian Htot read 

n 
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2a/i + + vTT^^ 



(|.9)i|5)2-|e)i|e) 
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1 + ^1 + '. 



2^1-^7^2 + ^1+7^2 

(|e)i|5>2-|5>i|e)2) 



(|e)i|5>2 + |g)i|e)2) 



(l5>i|.9>2 + |e>i|e)2) , 
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(|5)i|5>2-|e)i|e) 



2y^l+7/2 -^l+T^'^ 
1 _ /I + „2 

, ^ ^ (|e)i|5>2 + |g)i|e)2) 



2^14-7/2- yiT^ 



(4) 



with eigenvalues Ei = -2\fo?^~lP , E2 = -2 J, £'3 = 2J 



and Ei = 2Va2 + J2. 

It is reasonable to consider as initial state of the two 
atoms the ground state |.g)i|5)2; then we can expand it 
over the eigenstates basis (0J as 



\^{Q))^\g),\g)2^Y.^,\i^,), 



(5) 



with 



(1 - yiT^) ^i + 7y2 + ^TT^ 



277^1 + 7/2 



C2 = 0, 



C3 



1 

71' 



(1 + ^/l + W) ^l + r,-^ - ^/TTW 



277^1 + 7/2 



(6) 
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3 System dynamics = -i[iJtot,p] +D[c+]p + D[c_]p, (9) 



The evolution of the state © under Htot gives 



(7) 



where we have introduced the scaled time r = Jt. 

In Ref.^ it was shown that homodyne measurement 
of the light scattered by an atom allows indirect measure- 
ment of its spin flip operators. Then, let us consider, such 
type of local measurements so that after combining ho- 
modyne currents, the total current I{t) carries out infor- 
mation about the observable O = ctJ^'' — (T2^' . Its variance 
over the state JT)! is 



1 — cos 



(8) 



1 

Notice that this quantity being strictly less than 2 at al- 
most any time, shows the presence of correlations for the 
state JT)). On the other hand, in Ref. 5 it has been shown 
that the state lO exhibits entanglement at almost any 
time. We are thus led to ascribe the correlations of Eq.® 
to the presence of entanglement in Eq.(|7|, though this 
would not generally true. Then, we are going to consider 
the quantity O as a "marker" of entanglement while char- 
acterizing the open system dynamics. 

When we include the effect of spontaneous atomic de- 
cay at rate 7, the dynamics of the two distant atoms is 
described by the master equation 

p = -i [Htot.p] + D [cTi] p + D [(T2] p 



where c± = (cri ± 02)!^/^ and D is the Lindblad decoher- 
ence superoperator, i.e. D[a]& = aha^ — ah/2 — ha^a/2. 
The following replacements J/7 —^J,a/"f—^a,jt—ft 
have been made deriving Eq.Q. 

The steady state solution of Eq. can be easily found 
by writing the density operator in a matrix form, in the 
basis {|e)i|e)2, |ff)i|e)2, |e)i|5)2, l5>i|5>2}, as 



^ A Bi + iB2 Ci + iC2 Vi + W2 ^ 



Pss = 



(10) 



Bi — iB2 £ + iT2 Gi + iG2 

Ci - 162 Tx - iT2 H Ti+ 1X2 
- Qi - i'D2 1i - il2 C 
while the matrix representation of the other operators (in 
the same basis) comes from 



CTl 



'^O 0^ 



10 

10 



0-2 



^0 0^ 




10 
10 



(11) 



By inserting these matrices in the r.h.s. of Eq.© and 
equating to at l.h.s. we are left with a set of 16 linear 
equations from which we can calculate (together with the 
normalization condition tv{p) = 1) all the real coefficients 
of the matrix IjlOl) . namely 

1 



1 ^ 
Ci = -^8a^ 

1 



^2=0, 
C2=0, 

1 



:16aV, 



1 



£ = — (lea" -t- 
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-8aJ, 



with 



£ = — (16a'* + Sa^ + 1 + 16 J^) , 



S = 64a^ + IGa^ + 1 + 16 . 



(12) 



(13) 



4 Stationary entanglement 

One can use the concurrence as measure of the degree 
of entanglement between two qubit described by density 
operator p [5]. It is defined as 



C(p) = max {0,6-6-6-^4} 



(14) 



where 6's are, in decreasing order, the nonnegative square 
roots of the moduli of the eigenvalues of the non-hermitian 
matrix pp. Here p is the matrix given by 



(15) 



where p* denotes the complex conjugate. 

The stationary state concurrence C'o = C {pss ) is shown 
in Fig[21 It is clear that a relevant amount of entanglement 
persists at steady state only for large values of the coupling 
constant, i.e. J ^ 1, (when the original J is much greater 
than 7). 

5 Feedback action 

We can now think to stabilize the entanglement, i.e. to 
prevent its degradation, by using a feedback action on 




Fig. 2. Concurrence Co of the steady state plotted versus the 
driving strength a and the coupling constant J. 

the driving fields (LI and L2 of Fig^J accordingly to the 
measured quantity O which should reveal the status of 
nonclassical correlations. Then we act on the system with 
a local feedback operator 



F : 



^ ( iv) 



V2 



(16) 



where A represent the feedback strength (already scaled 
by 7, i.e. A/ ^ A). The choice of F is motivated by the 
fact that feedback mediated by indirect (homodyne) mea- 
surement requires, to squeeze the variance of a variable 
(O), a driving action on the conjugate variable [H]- 
The master equation then becomes [5] 



P = -i {lltot,p\ + D [c+] p 

+D[c_-iF]p-^ cLf^Fc^.P 



(17) 



In the above equation, the feedback operator F appears 
in the Hamiltonian term describing the driving effect, as 
well as inside the decoherence superoperator accounting 
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Fig. 3. Concurrence Cft of the steady state plotted versus 
the driving strength a and the coupling constant J in presence 
of feedback action. For each value of a and J, the feedback 
strength is chosen to be the optimal. 

for quantum noise carried back into the system from mea- 
surement. 

The master equation p7|) can be solved at steady state 
with the same method of Eq.lO, obtaining p(^. However, 
the analytical expression is quite cumbersome, hence not 
reported at all. The state p(g allows us to (numerically) 
calculate its concurrence. In particular, we have evaluate 
the quantity 



(18) 



that is shown in Fig|31 We can see that feedback improves 
the available entanglement with respect to previous case 
(Fig[2Jl. Feedback seems especially powerful at small val- 
ues of J where entanglement was very fragile (it somehow 
enforces the coupling effect). 



Fig. 4. Concurrences difference Cft — Co plotted versus the 
driving strength a and the coupling constant J. 

To better compare the results with and without feed- 
back, in Fig^we have shown the difference Cft — Cq. 



6 Conclusion 

We have shown the possibility to improve the steady state 
entanglement in an open quantum system by using a feed- 
back action. Although the improvement is not very high 
the above result represents a proof of principle about the 
possibility of controlling entanglement through feedback. 
A complementary possibility to increase entanglement be- 
tween atoms subject to joint measurements with feedback 
has been then proposed |lUj . 

Since our method only relies on Local Operations and 
Classical Communication (LOCC), what we have obtained 
is perhaps related to entanglement purification . 
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To improve the presented model one should find the 6. D. Gunlycke, S. Bose, V. M. Kendon and V. Vedral, Phys. 
best entanglement witness [T^] to measure, and then opti- Rev. A 64, (2001) 042302. 

mize the feedback action (operator). This can be phrased 7. H. M. Wiseman and G. J. Milburn, Phys. Rev. A 49, (1994) 
in terms of a numerical optimization problem and is left 1350. 

for future work. Moreover, since entanglement is a system W- K. Wootters, Phys. Rev. Lett 80, 2245 (1998). 

, , 1- -i ii r 11 1 1 1-1 i i 9. S. Mancini, D. Vitali and P. Tombesi, Phys. Rev. A 61, 

state peculiarity, other leedback procedures, like state es- i j > 

timation based feedback could be more powerful. 

Summarizing, although we have proved the possibility 
of feedback control of entanglement, its effectiveness re- 
mains difficult to quantify in nonlinear systems (like that 
studied). Probably, investigations in linear systems would 
be more fruitful. 
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